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ABSTRACT 

Massive merging black holes will be the primary sources of powerful gravi- 
tational waves at low frequency, and will permit to test general relativity with 
candidate galaxies close to a binary black hole merger. In this paper we identify 
the typical mass ratio of the two black holes but then show that the distance 
where gravitational radiation becomes the dominant dissipative effect (over dy- 
namical friction) does not depend on the mass ratio; however the dynamical 
evolution in the gravitational wave emission regime does. For the typical range 
of mass ratios the final stage of the merger is preceded by a rapid precession and 
a subsequent spin-flip of the main black hole. This already occurs in the inspiral 
phase, therefore can be described analytically by post-Newtonian techniques. We 
then identify the radio galaxies with a superdisk as those in which the rapidly 
precessing jet produces effectively a powerful wind, entraining the environmental 
gas to produce the appearance of a thick disk. These specific galaxies are thus 
candidates for a merger of two black holes to happen in the astronomically near 
future. 

Subject headings: compact binaries, gravitational radiation, radio galaxiesjets 
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1. Introduction 

The most energetic phenomenon that involves general relativity in the observable uni- 
verse is the merger of two supermassive black holes (SMBHs). Therefore the study of these 
mergers may provide one of the most stringent tests of general relativity even before the dis- 
covery and precise measurement of the corresponding gravitational waves (see, e.g., Schafer 
2005). 

Most galaxies have a central massive black hole (Kormendy & Richstone 1995, Sanders 
& Mirabel 1996, Faber et al. 1997), and after their initial growth (for one possible example 
how this might happen, see Munyaneza & Biermann, 2005, 2006), their evolution is governed 
by mergers. Therefore, the two central black holes also merge (see Zier & Biermann 2001, 
2002; Biermann et al. 2000; Merritt & Ekers 2002; Merritt 2003; Gopal-Krishna et al. 2003, 
2004, 2006; Gopal-Krishna & Wiita 2000, 2006; Zier 2005, 2006, 2007). Before the two black 
holes get close, the galaxies begin to round each other, distorting the shape of a radio galaxy 
fed by one or both of the two black holes; thence the Z-shaped radio galaxies (Gopal-Krishna 
et al. 2003). When they merge, under specific circumstances to be clarified in this paper, 
a spin-flip may occur. For a black hole nurturing activity around it, the spin axis defines 
the axis of a relativistic jet, and therefore a spin-flip results in a new jet direction: thence 
the X-shaped radio galaxies (Rottmann 2001, Chirvasa 2001, Biermann et al. 2000, Merritt 
and Ekers 2002). In fact, observations suggest that all activity around a black hole may 
result in a relativistic jet even for radio-weak quasar activity (Falcke et al. 1996, Chini et al. 
1989a, b). A famous color picture showing the past spin-flip of the M87 black hole (Owen 
et al. 2000) clearly shows a weak radio counter-jet, misaligned with the modern active jet 
by about 30°. The feature of the X-shaped radio galaxy jets is so common and yet very 
short-lived that all radio galaxies may have been through this merger (Rottmann 2001), and 
thus should have undergone a spin-flip. This can be also deduced from the observation that 
many compact steep spectrum sources show a misaligned double radio structure, where an 
inner pair of hot spots is misaligned with an outer pair of hot spots (Marecki et al. 2003). 
We conclude that theoretical arguments and observations consistently suggest that black 
holes merge and result in a spin-flip. 

From these and some other data we deduce a few basic tenets that the theory needs to 
explain: 

1. In the X-shaped radio galaxies the angles between two pairs of jets in projection are 
typically less than 30 degrees. The real angles can be even about 45°. The jets are believed 
to signify the spin axis of the more active (therefore presumably the more massive) black 
hole before the merger and the spin axis of the merged black hole. Therefore, a substantial 
spin-flip should have occurred. 
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2. In the X-shaped radio galaxies one pair of jets has a steep radio spectrum. This 
implies that it has not recently been resupplied energetically, it is an old pair of jets; and 
its synchrotron age is typically a few 10 7 years. The other pair of jets has a relatively flat 
radio spectrum (this is the new jet; Rottmann 2001). Radio continuum spectroscopy thus 
supports the spin-flip model. 

3. Again, as Rottmann (2001) shows, the statistics of X-shaped radio galaxies are 
such, that every radio galaxy may have passed through this stage during its evolution. This 
matches with arguments based on far-infrared observations that central activity in galaxies 
such as starbursts and feeding of the activity of a central black hole, is often, maybe always, 
preceded by a merger of galaxies (Sanders & Mirabel 1996). 

4. There is another critical observation of the spectrum of radio galaxies. For many 
of them the radio spectrum has a low-frequency cutoff, suggesting a cutoff in the energy 
distribution of the electrons at approximately the pion mass (the electrons / positrons are 
decay products from pions, produced in hadronic collisions; Falcke et al. 1995, Biermann et 
al. 1995, Falcke and Biermann 1995a, 1995b, 1999, Gopal-Krishna et al. 2004). Hadronic 
collisions with ensuing pion production at the foot ring of the radio jet occur naturally and 
thermally in the case that the rotation parameter of the black hole is larger than 0.95, and 
if the foot ring is an advection-dominated accretion flow (ADAF) or radiatively inefficient 
accretion flow (RIAF; Donea & Biermann 1996, Mahadevan 1998, Gopal-Krishna et al. 
2004). If this is true for all radio galaxies, the spin of the black hole both before and after 
the spin-flip must be more than 95 % of the maximally allowed value. This is a major 
constraint on the process of the spin-flip. If we assume that this holds for all radio galaxies, 
then a fortiori it also holds for those which have just undergone a binary black hole merger, 
and so their spin ought to be high as well. 

5. When two black holes merge, the emission of strong gravitational waves is certain 
(Peters & Mathews 1963, Peters 1964, Thorne 1979). Compact binaries are driven by grav- 
itational radiation through a post-Newtonian (PN) regime (the inspiral), a plunge and a 
ring-down phase toward the final state. It is commonly believed than the spin-flip phe- 
nomenon is likely to be caused by the gravitational radiation escaping the merging system 
(Rottmann 2001, Biermann et al. 2000, Merritt and Ekers 2002). Recent numerical work 
on the final stages of the coalescence supports this (see Briigmann 2007; Campanelli et al. 
2007a, b; Gonzalez et al. 2007a, b). 

Therefore, it is mandatory to investigate what happens when the two black holes get 
close to each other, and this we propose to treat in this paper. We present here a model which 
allows to have a merger transition going from a high-spin stage to another high-spin stage, 
using mostly physical insight from outside of the innermost stable orbit (ISO). In contrast 
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with available numerical simulation, our method, limited to a certain typical range of mass 
ratios of the two black holes, has the advantage that the evolution of the compact binary 
can be treated in the framework of an analytical PN expansion with two small parameters. 

In Section 2, we review the current state of observations on the masses of supermassive 
galactic black holes, which roughly scale with the bulge masses of their host galaxies. The 
observations suggest that the most massive black holes have about 3 x 10 9 solar masses 
(Mq) and the most reliable determination of the low-mass central black hole (in our galaxy) 
is about 3 x 10 6 M & (Ghez et al. 2005, Schodel & Eckart 2005). There is some evidence for 
central massive black holes of slightly lower mass (Barth et al. 2005), but the error bars are 
very large. This implies that the maximum mass ratio is about 10 3 . We carefully analyze 
the statistics and argue in Section 2 that mass ratios in the range 3 : 1 to 30 : 1 cover most of 
the plausible range in mergers of galactic central black holes. Roughly speaking, this means 
that typically one mass is dominant by a factor of order 10. Therefore, we find that neither 
the much discussed case of equal masses nor that of the extreme mass ratios (test particles 
falling into a black hole) describes typical central galactic SMBH mergers. 

In Section 3, we study the relative magnitudes of the spin of the dominant black hole 
and of the orbital angular momentum of the system. Their ratio depends on two factors: 
the mass ratio and the separation of the binary components (the inverse of which scales with 
the post-Newtonian parameter). We show that for the typical mass ratio interval the orbital 
momentum left when the system is reaching ISO is much smaller than the dominant spin. 
So in the typical mass range case whatever happens during the plunge and ring-down phases 
of the merger, in which the remaining orbital momentum is dissipated, it cannot change 
essentially the direction of the spin. By contrast, for equal mass mergers the orbital angular 
momentum dominates until the end of the inspiral, while for extreme mass ratio mergers the 
larger spin dominates from the beginning of the gravitational wave driven merger phase. 

In Section 4, we discuss the transition from the dynamical friction dominated regime to 
the gravitational radiation dominated regime, in order to establish the initial data for the PN 
treatment. The interaction of the black holes with the already merged stellar environment 
generates a dynamical friction when the separation of the black holes is between a few parsecs 
(pc) and one hundredth of a pc. Gravitational radiation has a small effect in this regime. 
Due to the dynamical friction, some of the orbital angular momentum of the binary black 
hole system is transferred to the stellar environment, such that the stellar population at the 
poles of the system tends to be ejected and a torus is formed (Zier & Biermann 2001, 2002, 
Zier 2006). This connects to the ubiquitous torus around active galactic nuclei (AGNs), 
detected first in X-ray absorption (Lawrence & Elvis 1982, Mushotzky 1982), and later 
confirmed by optical polarization of emission lines (Antonucci & Miller 1985). Dynamical 
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friction is enhanced as in a merger the phase-space distribution is strongly disturbed by large 
fluctuations of the mass distribution (Lynden-Bell 1967, Toomre & Toomre 1972, Barnes & 
Hernquist 1992, Barnes 2001). There had been a major worry that the two black holes 
stall in their approach to each other (Valtonen 1996, Yu 2003, Merritt 2005, Milosavljevic 
& Merritt 2003a, b, Makino & Funato 2004, Berczik et al. 2005, 2006, Matsubashi et al. 
2007) before they get to the emission of gravitational waves; that the loss-cone mechanism 
for feeding stars into orbits that intersect the binary black holes is too slow. However, Zier 

(2006) has demonstrated that direct interaction with the surrounding stars slightly further 
outside speeds up the process, and so very likely no stalling occurs. Relaxation processes 
due to cloud/star-star interactions are rather strong, as shown by Alexander (2007), using 
observations of our galaxy. These interactions repopulate the stellar orbits in the center of 
the galaxy. New work by Merritt et al. (2007) is consistent with Zier (2006) and Alexander 

(2007) . Also in a series of papers Sesana, Haardt, and Madau have recently shown that 
even in the absence of two-body relaxation or gas dynamical processes, unequal mass and/or 
eccentric binaries with the mass larger than 10 5 M & can shrink to the gravitational wave 
emission regime in less than a Hubble time due to the binary orbital decay by three-body 
interactions in the gravitationally bound stellar cusps (Sesana et al. 2006, 2007a, 2007b). 
Finally, Hayasaki (2008) has considered the "last parsec problem" under the assumption of 
the existence of three accretion disks: one around each black hole and a third one, which is 
circumbinary. The circumbinary disk removes orbital angular momentum from the binary 
via the binary-disk resonant interaction, however, the mass transfer to each individual black 
hole adds orbital angular momentum to the binary. The critical parameter of the mass 
transfer rate is such that for SMBH binaries, it becomes larger than the Eddington limit, thus 
these binaries will merge within a Hubble time by this mechanism. The angular momentum 
transfer from orbit to disk was already considered as a key physical concept in binary stars by 
Biermann & Hall (1973). All these recent works suggest that by one mechanism or another 
the SMBHs will approach each other to distances smaller than approximately one hundredth 
pc, when the gravitational radiation becomes the dominant dissipative effect. In Section 4, 
we analyze the characteristic timescales of the dynamical friction and gravitational radiation 
as function of the total mass, stellar distribution radius and mass ratio of the compact 
binary and we establish the values of the transition radius and PN parameter, for which the 
gravitational radiation is overtaking dynamical friction. 

In Section 5, we discuss the post-Newtonian evolution of the compact binaries, following 
Apostolatos et al. (1994) and Kidder (1995). The new element is the emphasis of the role of 
the mass ratio as a second small parameter in the formalism. The leading order conservative 
effect contributing to the change in the orientation of spins is the spin-orbit (SO) coupling. 
The backreaction of the gravitational radiation, which is the leading order dissipative effect 
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below the transition radius, appears at one PN order higher. We show here that for the 
characteristic range of mass ratios the spin-flip occurs during the gravitational radiation 
dominated inspiral regime, outside ISO. In the process we evaluate the timescales for the 
change of the spin tilt as compared to the timescales of precessional motion and gravitational 
radiation driven inspiral. As a by-product, we are able to show that for the typical mass 
range the so-called transitional precession occurs quite rarely. 

We interpret and discuss the resulting model in Section 6. Here, we give a tentative 
outline of the time sequence of the activity of two merging galaxies, leading to an AGN 
episode of the primary black hole. A recent review of the generic aspects of these galaxy 
nuclei as sources for ultra high energy cosmic rays is in (Biermann et al. 2008). 

Finally, we summarize our findings in the concluding remarks. Following our arguments 
about the phase just barely before the merger we propose there that the superwinds in radio 
galaxies (Gopal-Krishna et al. 2007, Gopal-Krishna & Wiita 2006) are in this stage, as the 
rapidly precessing jet acts just like a powerful wind. The primary goal of our paper is to put 
the derived physics into observational context, so as to allow tests to be done in radio and 
other wavelengths. 

2. The relevant mass ratio range 

In Lauer et al. (2007), the mass distribution of galactic central black holes is described, 
confirming earlier work, and also consistent with a local analysis (Roman & Biermann 2006). 
Arguments based on Haring & Rix (2004), Gott & Turner (1977), Hickson (1982), and Press 
& Schechter (1974) reasoning lead to a similar result, as does a recent observational survey 
(Ferrarese at al. 2006b). Wilson & Colbert (1995) also find a broken power law. The 
probability for a specific mass ratio is an integral over the black hole mass distribution, folded 
with the rate to actually merge (proportional to the capture cross-section and the relative 
velocity for two galaxies), e.g., isomorphic to the discussion in Silk & Takahashi (1979) for 
the merger of clumps of different masses. The black hole mass distribution &bh{Mbh), 
the number of massive central black holes in galaxies per unit volume, and black hole mass 
interval, can be described as a broken power law, from about m a ~3x 10 6 M to about 
mj,~3x 10 9 M , with a break near m+ ~ 10 8 M . The lower masses have been discussed 
in some detail by Barth et al. (2005). The values of m a , m b and imply that we have 
two mass ranges of a factor of 30 each. The masses above 10 8 M Q are rapidly becoming rare 
with higher mass, so that the lower mass range is statistically more important. That ratio 
range is then 1 : 1 to 30 : 1; while in the higher mass range the maximal range of the masses 
is also 30 : 1. 
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The mass of the central massive black hole scales with the mass of the spheroidal com- 
ponent, as with the total mass of a galaxy (the dark matter), see Benson et al. (2007). The 
rate of black hole mergers is some fraction of all mergers of massive galaxies. If, as argued 
by Zier (2006) the approach of the two black holes does not stall, then each merger of two 
massive galaxies will inevitably lead to the merger of the two central black holes. This is 
supported by the statistical arguments of Rottmann (2001), using radio observations, that 
all strong central activity in galaxies may involve a merger of two black holes. Therefore, 
observational evidence suggests that black holes do merge, and do so on the rather short 
timescales of AGNs. 

The interactions and mergers of galaxies clearly depend on the three angular momenta: 
the two intrinsic spins, and the relative orbital angular momentum, as well as on the initial 
distance and relative velocity of the two galaxies. Once all these parameters are given, the 
evolution is quite deterministic. The observations of Gilmore et al. (2007) strongly suggest, 
that the initial seed galaxies are today's dwarf elliptical galaxies, all of which are consistent 
with a lower bound to a common total mass of 5 x 10 7 M . This implies that all galaxies, 
and a fortiori all central black holes, have undergone very many mergers. 

The observations of Bouwens & Illingworth (2006) and lye et al. (2006) strongly suggest 
that much of this merger history happened earlier than redshift 6, perhaps mostly between 
redshifts 9 and 6. Each individual merger runs along a well-defined evolutionary track, but 
all of these mergers are completely uncorrelated with each other. Therefore, the ensemble of 
very many mergers can be treated statistically, and this is what we proceed to do, using the 
constant mass ratio between the spheroidal component of galaxies and their central black 
holes. We thus use the merger rate of galaxies as closely equivalent to the merger rate of the 
central black holes. 

The statistics of the mergers is given by the integral for the number of mergers N(q) per 
volume and time for a given mass ratio q, defined to be larger than unity. This merger rate 
is the product of the distribution of the first black hole with the distribution of the second 
black hole multiplied by a rate F. The latter in principle depends on both the cross-section 
and relative velocity of the two galaxies, the velocities however are not very different, as the 
universe is not old enough for mass segregation. The cross-section in turn depends on the 
two masses, thus F = F(q,m). If we integrate for all cases, in which the first black hole 
is less massive than the second black hole, we undercount by a factor of 2, and we have to 
correct for this factor. The general relationship is 




(1) 
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It is likely that the more massive black hole, and so the more massive host galaxy, will 
dominate the merger rate F, so that it can be approximated as a function of qm alone, and 
a power law behavior with F ~ q^ with £ > should be adequate for a first approximation. 
To estimate £ roughly we just observe, that dwarf spheroidals have a core radius of a few 
hundred pc (Gilmore et al. 2007), while our Galaxy has a core radius of about 3 kpc (Klypin 
et al. 2002), so a factor of 10 in radius (10 2 in cross-section) for a factor of about 10 4 in mass, 
thus the exponent is likely to be approximately 1/2; therefore a reasonable first estimate for 
any cross-section is £ = 1/2. In this instance we use the approximate equivalence of galaxy 
mergers with black hole mergers. 

As the black hole mass distribution has a break at = 30, we use (m) ~ m~ a for 
the first mass range, and $bh{™) ~ m _/3 for the second. For the range q from 1 to 30 we 
have as a dominant contribution 




The various models shown in Lauer et al. (2007) show that a range of values of a and j3 
is possible, with a ranging between approximately 1 and 2, and f3 from 3 to larger values. 
Benson et al.(2007) propose a ~ 0.65. We adopt here the approximate values for a and (3 
of 1 and 3, to be cautious, and for £ we adopt 1/2. With these values the above integrands 
are monotonically decreasing functions and the integrals are dominated by the lower limits. 
Thus, the four terms scale with q as q^~ a , q~ 1+a , q^~^, and again q^ 13 . 

Let us consider the four terms: the first term is small galaxies merging with small 
galaxies, and so not very interesting, as the cross-section is low. However, for this distribution 
the number of mergers in the mass ratio range 30 : 1 to 3 : 1 versus 3 : 1 to 1 : 1 is about 5. 
The more extreme mass ratios are more common. For the second term this ratio of mergers 
in the two mass ratio ranges is about 14. As this is massive galaxies merging with smaller 
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galaxies (above and below the break m*), this is the most interesting case, and also quite 
common. The third term is almost negligible, and the fourth term adds cases to the second 
term with more extreme mass ratios, beyond 30 : 1, and so emphasizes the large mass ratio 
range. 

So, among the relevant cases the rate of mergers of mass ratio of more than 3 : 1 to those 
with a smaller mass ratio is in the range of 5 : 1 to 14 : 1, about an order of magnitude. 
Focussing on those cases where one black hole is at 10 s M or larger, the ratio is larger 
than 14 : 1. Speculating that the exponent £ could be larger would enhance all these effects; 
enlarging a would weaken them. Therefore we will deal in the following with this much more 
common extended mass ratio range 30 : 1 to 3 : 1, which as will be shown, allows to use 
analytical methods. 



3. The spin and orbital angular momentum in the PN regime 

We assume the compact binary system to be composed of two masses m^, % — 1, 2 , each 
having the spin Si. By definition, the characteristic radius Ri of compact objects is of the 
same order of magnitude that the gravitational radius Rq = Grrii/c 2 (where c is the velocity 
of light and G is the gravitational constant). Therefore, the magnitude of the spin vector 
can be approximated ctS S{ 77ij RiVi Gm^Vi/c 2 , where Vi is the characteristic rotation 
velocity of the ith compact object As black holes rotate fast due to accretion, V*/c is of 
order unity. Equivalently we can introduce Si = (G/c)m 2 Xi, with Xi being the dimensionless 
spin parameter. Then maximal rotation implies Xi = 1- 

The PN expansion is done in terms of the small parameter 

Gm v 2 

(4) 



c z r cr 



where m = m 1 +m 2 is the total mass and v is the orbital velocity of the reduced mass particle 
\i = mim 2 /m, which is in orbit about the fixed mass m (according to the one-centre problem 
in celestial mechanics) . The two expressions for e are of the same order of magnitude due to 
the virial theorem. As in certain expressions odd powers otv/c may occur, it is common to 
have half-integer orders in the post-Newtonian treatment of the inspiral of a compact binary 
system. 

Whenever the masses of the two compact objects are comparable, either of Grrii/c 2 r 
also represent one post-Newtonian order. However, as we have argued before, for colliding 
galactic black holes it is typical that their masses differ by 1 order of magnitude, so that we 
have a second small parameter in the formalism. By choosing m 2 as the smaller mass, we 



-10- 



can also define the mass ratio 

i/ = ^ = -e(0,i) . (5) 

m 1 q 

In the literature the symmetric mass ratio 

is also frequently employed. The two mass ratios are related as 

and for small v we have rj = v — 2u 2 + O (u 3 ). 

For the typical mass ratio range of SMBH binaries either rj or v can be chosen as 
the second small parameter in the formalism. However, while these stay constant, the PN 
parameter e evolves during the inspiral toward higher values. Indeed, the separation of the 
components of the binary with m = 1O 8 M evolves as 

r=^ = £ = 4.7813x10"^, (8) 

c 2 e 2e e ' w 

where r$ represents the Schwarzschild radius. The interaction of the galactic black holes 
with the stellar environment begins when the black holes are a few kpc away from each 
other (then e ~ 10~ 8 ). The dynamical friction becomes subdominant at about 0.005 pc 
(Zier & Biermann 2001, Zier 2006), when the gravitational radiation becomes the leading 
dissipative effect. Thus, e = e* ~ 10~ 3 is the value of the PN parameter for which the 
gravitational radiation is driving the dissipation of energy and orbital angular momentum. 
Then follows the inspiral stage of the evolution of compact binaries, which continues until 
the domain of validity of the post-Newtonian approach is reached, at few gravitational radii, 
at ISO. Further away a numerical treatment is necessary in order to describe the plunge, 
which is finally followed by the ring-down. The PN formalism can be considered valid until 

e w nr 1 . 

Theoretically, it is possible for a small u, that at certain stage of the inspiral, the 
increasing e 1 ^ 2 becomes of the same order of magnitude as v and later on it even exceeds v. 
Such a situation would shift the numerical value of several contributions to various physical 
quantities into the range of higher or lower PN orders, depending of the involved power of 



The spin ratio (for similar rotation velocities V\ fa V 2 ) can be expressed as 

S2 ( m 2 \ 2 2 
Si \rrij 



(9) 
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Table 1: The evolution of the ratio S±/L ~ e x l 2 v 1 in the range e — 10 3 ^-10 
for various values of the mass ratio v. 







e « 10~ 3 




e « 10" 1 


v=\ 


0.03 


(Si « L) 


0.3 


(^i < L) 


v= 1/3 


0.1 


(Si < L) 


1 


(ft « L) 


v = 1/30 


1 


(Si « L) 


10 


(5i > L) 


v = 1/900 


30 


(Si » L) 


300 


(Si » L) 



The ratio of the spins to the orbital angular momentum becomes 

te^v, (10) 



S 2 Gm 2 V 2 /c 2 f Gm\ ( c\ fV 2 \ m 2 y 2 



L jirv \ c 2 r J \v/ \ c J m x 

e 1 ' 2 ^ 1 . (11) 



Si S 2 Si 



L L S- 



2 



We note that the approximations in the above formulae (I9l- ([TT1) are related only to the fact 
that we have assumed maximal rotation (thus Vi/c < 1). First, we note that the above ratios 
involving the spins of the compact objects already contain e 1 ^ 2 . Thus, the counting of the 
inverse powers of c 2 is not equivalent with the PN order, when compact objects are involved. 
Further, while the ratio S 2 /L is shifted toward higher orders by a small v (therefore S 2 ^ L 
during all stages of the inspiral), the order of the ratio of the spin of the dominant black hole 
to the magnitude of the orbital angular momentum is not fixed. Indeed, it is determined by 
the relative magnitude of the small parameters e and v. As e increases during the inspiral, 
whenever v falls in the range of e 1//2 , the initial epoch with S\ < L is followed by S\ ~ L 
and Si > L epochs (Tabled]). 

We have concluded in the previous section that the range of mass ratios q between 
3 : 1 and 30 : 1 is the most common. For such binaries the sequence of the three epochs 
Si < L, Si ~ L, and S± > L is fairly representative. We call this intermediate mass ratio 
mergers, which has to be contrasted with the case of equal mass mergers, where the orbital 
angular momentum dominates throughout the inspiral; and with the case of extreme mass 
ratio mergers (which we define as having mass ratios larger than 30 : 1), where the larger 
spin dominates from the beginning of the inspiral to the end of the PN phase, as can be seen 
from our Table [TJ 
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4. The timescales 

The value e ~ 10~ 3 from which the PN analysis with the gravitational radiation as the 
leading dissipative effect can be applied was adopted in the previous section for a compact 
binary with total mass m = 10 8 M G and mass ratio v = 10 _1 . This was based on the analysis 
in Zier & Biermann (2001) and Zier (2006), where it was shown that at around 5 x 10~ 3 
pc gravitational radiation takes over from dynamical friction in the interaction with stars in 
the angular momentum loss of the black hole binary. Further arguments for the binary to 
reach the gravitational wave emission regime were presented by Alexander (2007), Sesana et 
al. (2006, 2007a, 2007b), and Hayasaki (2008). 

In this section we raise the question, whether the value of the transition radius (and the 
corresponding value of the PN parameter) depends on m and v. In order to answer this, we 
compare the characteristic timescales of the gravitational radiation and dynamical friction. 

The timescale of gravitational radiation (as will be derived in Section 5) is 

1 L 32c 3 _ 4 



t gw L 5Gm 



e* V . (12) 



The timescale for the secondary black hole to lose angular momentum by gravitational 
interaction with the surrounding stellar distribution is (Binney & Tremaine 1987) 

v 3 fAv\ 2 . . 

tfr = 7^— T — • (13) 



2<KG 2 m 2 Pdistr^ V v J 

With the maximally allowable change in the velocity Av/v = 1 we find the relevant full 
timescale. Here, A = ln(b max /bmin) is the logarithm of the ratio of the maximal distance 
within the system, divided by the typical distance between objects. The latter is large for 
clouds, so the ratio is low and A of the order unity, while for stars or dark matter particles A 
can be taken as 10 20. For the merger, the estimate based on clouds is more appropriate, 
thus, following the reasoning of Binney & Tremaine (1987) we adopt A = 3. The compact 
stellar distribution with density pdistr, radius rdistr, and mass nidistr is of the same order in 
mass as the black hole binary (Zier & Biermann 2001, Ferrarese et al. 2006a) with a scale 
r distr of a few pc, and so under the assumption of a spherically symmetric distribution we 
set 

^ i 1 A\ 

Pdistr = "j 3 • 14 

By employing the definitions of the PN parameter and mass ratio, we get 

1 9G 2 m 2 



tf r 2c 3 r 



, e- 3/2 r] (1 + u) . (15) 



distr 
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This gives the full timescale of the dynamical friction. 

The two timescales become comparable for a PN parameter: 



Gm \ 6/n 

K{y)[^-) , (16) 



v C T distr , 

with K (is) being a factor of order unity, defined as: 

2/11 

G (0.938, 1.064) . (17) 



64 (1 + Z/) 



K(u) = 

corresponding to the distance r* 

r . = A - lM (^) 5/1I r «m. (18) 

Notably, the dependence on the mass ratio is rather weak and in practice it can be neglected. 
Inserting then for r distr = 5 pc and using as a reference value for the mass m = 10 8 M Q , we 
obtain e* ~ 10~ 3 and r* ps 0.005 pc in agreement with the discussion in Zier & Biermann 
(2001). We also note that in fact any other reasonable value for A and Av/v will give a 
factor of order unity in Eq. (14), as this number arises by taking the power 2/11. The weak 
dependence on v through K is due to the same reason. 

The value of the PN parameter thus scales with m 6//u and radius r^ 6 /^ 1 as 

^ « io- s f-^rV /U f— V /U > (19) 

while the transition radius scales with m 5 / 11 and r G J^ r as 

/ m \ 5/n / r^str \ 6/1 1 

* V1° 8M ©/ V 5 P c 
For m = 1O 9 M and for the same radius of stellar distribution then r* « 0.01 pc. 

We conclude that both the dependence of the transition radius and the corresponding 
PN parameter on the total mass and on the stellar distribution radius are weak, while there 
is practically no dependence on the mass ratio. 



0.005 pc — — ^= . (20) 

y \10 8 M Q ) \5 pc) K ' 



5. The inspiral of spinning compact binaries in the gravitational radiation 

dominated regime 

In the first subsection of this section we present the conservative dynamics of an isolated 
compact binary in a post-Newtonian treatment, emphasizing the role of the second small 



-14- 



parameter, as a new element. Then in the second subsection we take into account the effect 
of gravitational radiation, deriving how the spin-flip occurs for the typical mass ratio range. 
The limits of validity of our results obtained by using these two small parameters will be 
considered below in subsection 5.3. 



5.1. Conservative dynamics below the transition radius 

The interchange in the dominance of either L or Si has a drastic consequence on the 
dynamics of the compact binary. To see this, let us summarize first the conservative dynam- 
ics, valid up to the second post-Newtonian order. The constants of the motion are the total 
energy E and the total angular momentum vector J = L + Si + S2 (Kidder et al. 1993). 
The angular momenta L, S ; are not conserved separately. The spins obey a precessional 
motion (Barker & O'Connell 1975, Barker & O'Connell 1979): 

S; = ft-, x Si , (21) 

with the angular velocities given as a sum of the spin-orbit, spin-spin, and quadrupole- 
monopole contributions. The latter come from regarding one of the binary components as a 
mass monopole moving in the quadrupolar field of the other component. 

The leading order contribution due to the SO interaction (discussed in Kidder et al. 
1993, Apostolatos et al. 1994, Kidder 1995, Ryan 1996, Rieth & Schafer 1997, Gergely at 
al. 1998a, 1998b, 1998c, O'Connell 2004), cause the spin axes to tumble and precess. The 
spin-spin (Kidder 1995, Apostolatos 1995, Apostolatos 1996, Gergely 2000a, 2000b), mass 
quadrupolar (Poisson 1998, Gergely & Keresztes 2003, Flanagan & Hinderer 2007, Racine 
2008), magnetic dipolar (Ioka & Taniguchi 2000, Vasuth et al. 2003), self-spin (Mikoczi et 
al. 2005) and higher order spin-orbit effects (Faye et al. 2006, Blanchet et al. 2006) slightly 
modulate this process. 

The SO precession occurs with the angular velocities 

n. ^ . (») 

where Ln = fir x v is the Newtonian part of the orbital angular momentum. The total 
orbital angular momentum L also contains a contribution Lso (Kidder 1995), which for 
compact binaries is of the order of e 3 / 2 Ln ■ 



Due to the conservation of J, the orbital angular momentum evolves as 
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° [(4 + 3u) Si + (4 + 3^) S a ] xL . (24) 



2c 2 r 3 

(By adding a correction term of order e 3 ^ 2 relative to the leading order terms, we have 
changed L N into L on the right-hand side of the above equation.) 



To leading order in v we obtain: 



1C 

Sx = -— L x Si , (25) 



2C 

L = ^S lX L. (26) 

(Again, an Lso term was added to Ln on the right-hand side of Eq. ff25l) . in order to have 
a pure precession of Si.) 

Thus, the leading order conservative dynamics gives the following picture: the dominant 
spin Si undergoes a pure precession about L, while L does the same about Si. Despite the 
precession ( 1231) . the spin S 2 can be ignored to leading order, as its magnitude is v 2 times 
smaller than Si, e.g., Eq. ([9]). By adding the vanishing terms (2G/c 2 r 3 )Si x Si and 
(2G/c 2 r 3 ) L x L to the right-hand sides of Eqs. (1251) and fT2"6"j) . respectively, we obtain 

2C 

Si = ^JxSi, (27) 
2C 

L = -JxL. (28) 

Thus, the precessions can also be imagined to happen about J, which represents an invariant 
direction in the conservative dynamics up to 2PN. 

Higher order contributions to the conservative dynamics slightly modulate this pre- 
cessional motion. In fact, for both the spin-spin and quadrupole-monopole perturbations 
an angular average L can be introduced, which is conserved up to the 2PN order (Gergely 
2000a). As L differs from L just by terms of order 2PN, and L is conserved, the real evolution 
of L differs from a pure precession only slightly. 

Finally, we note that as the SO precessions are 1.5 PN effects and the gravitational ra- 
diation appears at 2.5 PN, at the transition radius the SO precession timescale is e _1 times 
shorter than the timescale of dynamical friction. The modifications induced by the preces- 
sions in the angular momentum transfer toward the stellar environment will be discussed 
elsewhere (Zier et al. 2009, in preparation). 
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Fig. 1. — The old jet points in the direction of the original spin Si. When the two black holes 
approach each other due to the motion of their host galaxies, a slow precessional motion of 
both the spin and of the orbital angular momentum L begins (left figure) about the direction 
of the total angular momentum J, which is due to the spin-orbit interaction. Gravitational 
radiation carries away both energy and angular momentum from the system, such that the 
direction of J stays unchanged. As a consequence the precessional orbit slowly shrinks and 
the magnitude of L decreases. This is accompanied by a continuous increase in the angle a 
and a decrease in f3. When the magnitudes of L and S become comparable (middle figure), 
the precessional motions are much faster (for typical values see Table [2]). In the typical 
mass ratio range v = 1/3 -f- 1/30 the magnitude of L becomes small as compared to the 
magnitude of the spin, which is unchanged by gravitational radiation (except the upper 
boundary of the range at v — 1/3 when L and S are still comparable). Before reaching the 
innermost stable orbit, the spin becomes almost aligned to the (original direction of the) 
total angular momentum, and a new jet can form along this direction. Therefore, for the 
typical mass ratio range the spin-flip phenomenon has occurred in the inspiral phase and 
not much orbital angular momentum is left to modify the direction of the spin during the 
plunge and ring-down. In the regime in between the initial and final states the precessing 
jet acts as a superwind, sweeping away the environment of the jets. 
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5.2. Dissipative dynamics below the transition radius 

Dynamics becomes dissipative from 2.5 PN orders. Then gravitational quadrupolar ra- 
diation carries away both energy and angular momentum. Orbital eccentricity is dissipated 
faster, than the rate of orbital inspiral (Peters 1964), thus the orbit will circularize^] Con- 
sidering circular orbits and averaging over one orbit gives the following dissipative change in 

L- = -^f^V /2 L, (29) 
5r \ err / 

where L represents the direction of the orbital angular momentum. Then the total change 
in L is given by the sum of Eqs. ( 1261) and ( 1291) . 



The spin-induced quasi-precessional effects both modulate the dynamics and they have 
an important effect on gravitational wave detection (see Lang & Hughes 2006, 2008, Racine 
2008, Gergely & Mikoczi 2008). 

The dissipative dynamics, with the inclusion of the leading order SO precessions and the 
dissipative effects due to gravitational radiation, averaged over circular orbits, was discussed 
in detail in Apostolatos et al. (1994) for the one-spin case S2 = and equal masses v — 1. 
For the typical mass ratio v e (1/30, 1/3), and keeping only the leading order contributions 
in the //-expansion also gives S2 = (the leading order contributions in S2 are of order 
z/ 2 ). In this subsection we will analyze in depth the angular evolutions and the timescales 
involved. 



As for any vector X with magnitude X and direction X one has X = XX + XX, the 
change in the direction can be expressed as X = ^X — Xx] / X. Also the identity X 2 = X 2 

gives X = X • X. Then Eqs. (EZD"® imply 

Si = 0, 



1C 

Si = -^JxSx 



c 2 r 3 



L 



32G> 2 (Gm\ 
5r V c 2 r / 



5/2 



1C 

L = -JxL. (30) 



c 2 r 3 



*As we focus here on spin-flips, we will not dwell on possible recoil as a result of the momentum loss due 
to gravitational radiation in the merger of the two black holes (see for example Brugmann 2008, Gonzalez et 
al. 2007a, b). The accuracy of determining the distance between two separate and independent active black 
holes (Marcaide & Shapiro 1983, Brunthaler et al. 2005) is reaching a precision, which may soon allow for 
recoil to be measurable; no such evidence has been detected yet. 
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The total angular momentum J is also changed by the emitted gravitational radiation. As 
no other change occurs up the 2PN orders, J =LL and 



J = HL-J 

L 



3 j 



L- L J J 



(31) 



Note that from the second Eq. fl3T|) it is immediate that the direction of J changes violently, 
whenever J is small compared to L. 

To leading order in v , the vectors L, Si, and J form a parallelogram, characterized by 
the angles a = cos -1 ■ and (3 = cos -1 ^Si-jj. From Eqs. fl30|) and fl3T|) . we obtain 

L . . 

a = sin«>0, (32) 

J 

$ = - sin a < . (33) 

In the latter equation, we have used that Si-L = cos [a + (3). 

Thus, we have found the following picture for the inspiral of the compact binary after the 
transition radius. By disregarding gravitational radiation, the SO precessions fl25|) and (1261) 
assure that the vectors L and Si are precessing about J (a fixed direction), but also about 
each other (then the respective axes of precession evolve in time). Gravitational radiation 
slightly perturbs this picture. The angle a + (3 between the orbital angular momentum and 
the dominant spin stays constant during the inspiral, even with the gravitational radiation 
taken into account. By contrast, the angle between J and L continuously increases, while 
the angle between J and Si decreases with the same rate. This also means that due to 
gravitational radiation, the vectors L and Si do not precess about J any more in an exact 
sense. They keep precessing about each other, however. 

The change in the total angular momentum J =LL is about the orbital angular momen- 
tum, which in turn basically (disregarding gravitational radiation) undergoes a precessional 
motion about J. This shows that the averaged change in J is along J. This conclusion, 
however, depends strongly on whether the precessional angular frequency fl p is much higher 
than the change in the angles a and (3. Indeed, if these are comparable, the component 
perpendicular to J in the change J =LL will not average out during one precessional cycle, 
as due to the increase of a it can significantly differ at the beginning and at the end of the 
same precessional cycle, see FigdJ 

The regime with Q p ^> a can be well approximated by a precessional motion of both 
L and Si about a fixed J, with the magnitudes of L and J slowly shrinking, the angle 
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a slowly increasing and (3 slowly decreasing. As a result, during the inspiral, the orbital 
angular momentum slowly turns away from J, while Si slowly approaches the direction of 
J. This regime is characteristic for the majority of cases, and it was called simple precession 
in Apostolatos et al. (1994). 

Whenever Q p rj a, the conclusion of having J in the direction of J does not hold for the 
average over one precession. This results in a change in the direction of J in each precessional 
cycle. The evolution becomes much more complicated (in fact no approximate analytical 
solution is known), and it was called transitional precession in Apostolatos et al. (1994). 

Let us see now when the two types of evolution typically occur. For this, we note that 
the inspiral rate L/L is of the order 

L 32 a < i 

L " "5G^ £ " ■ (34) 

while the precessional angular velocity Q p = 2GJ / c 2 r 3 gives the estimate 

n p ^< L . (35) 

Finally, the tilt velocity of L is of the order 

a « 32c E 7/2^ u (A s i n ( a + pj 



5Gm \ J 

<3ZC 9/2„,,-l / ^ 



r)G ^ e / V ( j J sin(a + /5) (36) 



S L 
sin a = -j sin (a + fi) w e xl2 v~ x — sin (a + (3) (37) 



We have used 



for the first and second expressions of a , respectively, together with Eq. ([TT|). For compar- 
ison, these are all represented in Table [2] for all three regimes characterized by S\/L w 0.3, 
Si ~ L and Si/L ~ 3, respectively. 

The numbers from the second line of Table [2] demonstrate that for the chosen typical 
example the precession timescale can get as short as a day, going from 3000 years to three 
years to a day in the three columns above. This last stage is obviously quite close to the 
plunge. From the first line we can infer upper limits of how close the merger actually is, so 
30 million years in Column 1, 300 years in Column 2, and a few months in Column 3. As the 
inspiral rate increases in time, rather than being a constant, these numbers are higher than 
the real values. The accuracy of the third estimate is further obstructed by the fact that 
after ISO the plunge follows, but as this comprises only a few orbits, the PN prediction can 
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Table 2: Order of magnitude estimates for the inspiral rate L/L, angular pre- 
cessional velocity Q p and tilt velocity a of the vectors L and Si with respect to 
J, represented for the three regimes with L > Si, L » Si and L < Si, charac- 
teristic in the domain of mass ratios v = 0.3 -j- 0.03. The numbers in brackets 
represent inverse time scales in seconds -1 , calculated for the typical mass ratio 
v = 10 _1 , post-Newtonian parameter 10~ 3 , 10 -2 and 10 _1 , respectively and 
m = 1O 8 M (then c 3 /Gm = 2 x 10~ 3 s" 1 ). 



L > Si L « Si L < Si 

-L/L (-ic) hps («io- n ) m^i (-lo- 7 ) 



be considered relevant as an order of magnitude estimate. By multiplying the numbers in the 
third line with the precession timescales Q~ 1 we actually get the relevant tilt angle, varying 
from 2 arcsec during one precession (6 x 10~ 4 arcseconds per year) at large separations to 3 
arcmin per precession (per day) close to the ISO. 

We see, that the rate of precession and the tilt velocity become comparable in the Si ~ L 
epoch (in which e 1 I 2 ij~ 1 m 1) for 



16 . . 
— sin [a + p) 
5 



1 « (eVY /S « e 1/a « * , (38) 



that is, for the chosen numerical example v = 10 _1 and for the square bracket of order unity, 
this gives J/L ~ 10 _1 and the rate of a ~ fl p ~ 10 -9 (this is still 100 times faster than the 
rate of inspiral). 

The total angular momentum J can be that small only if L and Si are almost perfectly 
anti-aligned, thus a + [3 ~ n — 5, 5 1, and L ~ Si- What does the condition for 
transitional precession (1381) mean in terms of their angle, how close should that be to the 
perfect anti-alignment? To answer the question we note that 

J L cos a + Si cos /3 

— = ~ o sin a , (39) 

LJ Lj 

Then, comparing with the estimates fl37j) and (1381 we conclude, that transitional precession 
can occur only if the deviation from the perfect ant i- alignment is of the order of z/ 3//2 . This 
is a highly untypical case in galactic black hole binaries. 
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5.3. The limits of validity 

One might seriously question whether pushing the values of the parameters beyond the 
range for which we use them would actually demonstrate that our approximations cannot 
possibly be correct. We would like to address the following two concerns specifically. 

a) Is the high orbital angular momentum limit L » Si, obtained by a sufficient increase 
of the separation of the two black holes a correct limit? 

b) Is the extreme mass limit v — > 0, for which the treatment of Section 5 may seem 
increasingly accurate, correct? 

Concerning the limit (a): in any expansion with a small parameter one condition always 
holds: the parameter must be small, and in any PN expansion we can always reach a stage, 
for which the physics basis fails as inadequate, as other additional physical processes become 
dominant, or for which there is no observational support, despite the mathematical validity 
of the expansion. 

From among the effects affecting the directions of the spins the dynamics discussed in 
Section 5 takes into account (1) the leading order conservative effect, given by the precession 
due to spin-orbit coupling and (2) the leading order dissipative effect due to gravitational 
radiation. Other conservative and dissipative effects are neglected, being weaker. Meaningful 
results can be traced from this model only when these assumptions hold. This implies that 
the post-Newtonian parameter e varies between 10~ 3 and 10 _1 , corresponding to orbital 
separations of 500 Schwarzschild radii = 0.005 pc to 5 Schwarzschild radii, given a 10 8 M & 
black hole. This is emphasized in the paragraph following Eq. (jSJ). 

The "initial" and "final" phases in the dynamics described above therefore refer to a 
well-defined range of orbital separation, they are not arbitrary. The choice for the initial state 
is further justified by the discussion of Section 4. One cannot apply the dynamics discussed 
above at arbitrarily large distances, where the orbital momentum indeed would dominate, 
simply because the dynamics is no longer valid there. At larger separations the leading 
order dissipative effect is due to dynamical friction, thus the discussion of the previous two 
subsections does not apply. 

Concerning (b): according to the summary presented in Table 1, there are three possi- 
bilities in the PN regime, where the dynamics discussed above holds: (i) from mass ratios 
v from 1 to 1/3 the orbital angular momentum dominates throughout the whole range of 
separations between 0.005 pc and 5 Schwarzschild radii, as noted; (ii) from mass ratios v 
from 1/3 to 1/30 initially the orbital angular momentum dominates over the spin, but their 
ratio is reversed at the final separation; (iii) from mass ratios v smaller than 1/30 the spin 
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is dominant throughout the process. 

Our claim is that the spin-flip is produced only if the total angular momentum (whose 
direction stays unchanged) initially is dominated by the orbital angular momentum, finally 
by the spin, thus only for the case (ii). 

In the dynamics presented in the previous two subsections we have neglected the second 
spin. As even for the highest mass ratio v — 1/3 in the regime (ii) the second spin is 1 
order of magnitude smaller than the leading spin, we consider this assumption justified. 
With decreasing mass ratios it becomes increasingly accurate to neglect the second spin, as 
according to Eq. (7) the ratio of the spins goes with v 2 . 

However, not all results of the previous subsection become increasingly accurate with 
a decreasing mass ratio. We emphasize, what is different in case (iii) as compared to (ii). 
The difference is in the initial conditions, which allow to obtain a spin-flip in case (ii), but 
not in (iii). Mathematically, the difference between these two cases can be seen from Eq. 
(33), showing that the angular tilt velocity of the dominant spin scales with z/ 2 . For the 
extreme mass ratios v < 1/30, e.g., Table 1, the spin dominates over the orbital angular 
momentum throughout the whole PN regime. Therefore, the ratio Si/ J is of order unity. 
With decreasing v however the change in the direction of the spin, represented by a goes 
fast to zero, thus no spin-flip is produced in the PN regime for extreme mass ratios. 

At the end of this subsection we derive an analytical expression relating a to the con- 
served a + (3, the evolving PN parameter e and the mass ratio v. From the Figure [1] 
J = L cos a + Si cos (3. By introducing the angle a + (3 and employing the estimate (ITT]) we 
obtain 

-«[1 + e 1/2 v~ 1 cos (a + (3)] cos a + e l/2 v~ l sin (a + (3) sin a . (40) 

Inserting this into the second expression ( 13T1) and rearranging we find 

sin 2a ^ sin (a + 0) 
1 + cos 2a ~ e~V 2 u + cos (a + (3) ' 

For an initial configuration of 0.005 pc (such that e = e* = 10 -3 ) and mass ratio v = 10 _1 , 
the initial misalignment between L and J is ai n ui a i ~ 18°, 10°, 0° for the dominant spin in 
the plane of orbit, spanning 45° with the plane of orbit and perpendicular to the plane of 
orbit (such that a+(3 = 90°, 45°, 0°), respectively. Then (3 initial = 72°, 35°, 0°. For the same 
mass ratio and relative configurations, the angle a at the end of the PN epoch (at e = 10 _1 ) 
becomes afi na i ~ 73°, 35°, 0°, respectively. This can be translated into a misalignment 
between Si and J of /3fi na i = 17°, 10°, 0°, and a spin-flip of A/5 = 55°, 25°, 0°, respectively. 
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5.4. Summary 

In the typical range of mass ratios v = 0.03 -r 0.3 the initial condition L > Si is always 
transformed into Si < L, but the transition is very rarely accompanied by the so-called 
transitional precession. In all other cases the precession is simple. As the precession angle 
of the dominant spin is decreasing in time from the given initial value to a small value, the 
precessional cone becomes narrower in time. At the end of the inspiral phase the dominant 
spin Si will point roughly along J. This means that a spin-flip has occurred during the post- 
Newtonian evolution, already in the inspiral phase of the merger. On the other hand, as the 
inspiral phase ends with L < S\, irrespective of what happens in the next phases, during the 
plunge and ring-down, L is not high enough to cause additional significant spin-flip. 

For smaller mass ratios (for extreme mass ratio mergers) the larger spin already domi- 
nates the total angular momentum from the beginning of the inspiral, thus no spin-flip will 
occur by the mechanism presented here. Alternatively, from the second expression f[3"6"j) one 
can see that the rate of tilt of the spin decreases with z/ 2 , thus it goes fast to zero in the 
extreme mass ratio case. However, as we argued in Section 2, such mass ratios are less 
typical for galactic central SMBH mergers. This also shows that an infalling particle will 
not change the spin of the supermassive black hole. 

For the (again untypical) equal mass SMBH mergers the orbital angular momentum 
stays dominant until the end of the inspiral phase. In this case, however the possibility 
remains open to have a spin-flip later on, during the plunge. 

6. Discussion 

The considerations from this paper lead to the following time sequence for the transient 
feeding of a SMBH including a merger with another SMBH. 

First: Two galaxies with central black holes approach each other to within a distance 
where dynamical friction keeps them bound, spiraling into each other. If there is cool gas 
in either one, it can begin to form stars rapidly, along tidal arms. The galactic central 
supermassive SMBH binary influences gas dynamics and star formation activity also in the 
nuclear gas disk, due to various resonances between gas motion and SMBH binary motion 
(Matsui et al. 2006), creating some characteristic structures, such as filament structures, 
formation of gaseous spiral arms, and small gas disks around SMBHs. If either galaxy 
happens to have radio jets, then due to the orbital motion, these jets get distorted and form 
the Z-shape (Gopal-Krishna et al. 2003, Zier 2005). 
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Second: The central regions in each galaxy begin to act as one unit, in a sea of stars and 
dark matter of the other galaxy. During this phase, as the cool gas from the other partner 
typically has low angular momentum with respect to the receiving galaxy, the central region 
can be stirred up, and produce a nuclear starburst (Toomre and Toomre 1972). The central 
black hole can get started to be fed at a high rate, but its emission will be submerged in 
all the far-infrared emission from the gas and dust heated by the massive stars produced in 
the starburst. In this case, there is dynamical friction, which can act so as to select certain 
symmetries, such as corotation, counter rotation, or rotation at 90° (as in NGC2685, a polar 
ring galaxy; Richter et al. 1994). 

Third: The black holes begin to lose orbital angular momentum due to the interaction 
with the nearby stars (Zier and Biermann 2001, 2002). Other mechanisms for angular 
momentum loss are also known (Sesana et al. 2006, 2007a,b, Alexander 2007, Hayasaki 2008). 
The two black holes approach each other to that critical distance where the interaction with 
the stars and the gravitational radiation remove equivalent fractions of the orbital angular 
momentum. Then, as shown in this paper, the spin axes tumble and precess. This phase 
can be identified with the apparent superdisk, as the rapidly precessing jet produces the 
hydrodynamic equivalent of a powerful wind, by entraining the ambient hot gas, pushing the 
two radio lobes apart and so giving rise the a broad separation (Gopal-Krishna et al. 2003, 
2007, Gopal-Krishna & Wiita 2006). Gopal-Krishna & Wiita (2006) emphasize the apparent 
asymmetry, which we propose to attribute to line-of-sight effects and the distortion due to 
the recent merger. The base of the radio structure is so broad and so asymmetric, that the 
central AGN will appear to be offset from the projected center of gap. The recent arguments 
of Worrall et al. (2007) seem to be consistent with this point of view. The spin direction of 
the combined two black holes is preserved, although the strength of the spin decreases. As 
during simple precession the total angular momentum shrinks considerably, but its direction 
is conserved, on the other side the magnitude of the spin stays constant, this means that the 
orbital angular momentum shrinks. For comparable mass binaries it will be still higher than 
the spin at ISO (therefore the dynamics below ISO, which can be analyzed only numerically, 
should be responsible for any spin-flip in the comparable masses case). For extreme mass 
ratio binaries the result of the shrinkage of the orbital angular momentum is L < S at ISO. 
Therefore, the spin at ISO should be roughly aligned with the direction of J = L + S, which 
(as initially L was dominant), is close to the direction of the initial L. 

In certain cases, especially for equal masses of the two black holes, a strong recoil has 
been found (Gonzalez et al. 2007a, b). However, as we noted earlier, the equal mass ratio is 
untypical. 

Fourth: The two black holes actually merge, and the merged black hole keeps the 
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spin axis from the orbital angular momentum of the previously existing binary, whenever 
the mass ratio is relatively large. In the case that the mass ratio is between 1 : 1 and 
1:3, then even at the innermost stable orbit a substantial fraction of the orbital angular 
momentum can survive, possibly leading to a spin-flip later on. This very short phase should 
be accompanied by extreme emission of low-frequency gravitational waves. The final stage 
in this merger leads to a rapid increase in the frequency of the waves, called "chirping", 
but this chirping will depend on the angles involved. The angle between the orbital spin of 
the combined two black holes, and intrinsic spin of the more massive black hole influences 
the highest frequency of the chirp; for a large angle this frequency will be lower than for a 
small angle between the two spins. Whether there is another observable feature, such as the 
induced decay of heavy dark matter particles, from the merger of the two black holes at that 
event such as speculated by Biermann & Frampton (2006) is not clear at this time. 

Fifth: Now the newly oriented more massive merged black hole starts its accretion disk 
and jet anew, boring a new hole for the jets through its environment. This stage can be 
identified perhaps with giga hertz peaked radio sources (GPS). If the new jet points at the 
observer, then 3C147 may be one example (Junor et al. 1999). 

Sixth: The newly oriented jets begin to show up over some kpc, and this corresponds to 
the X-shaped radio galaxies, while the old jets are fading but still visible. This also explains 
many of the compact steep spectrum sources, with disjoint directions for the inner and outer 
jets. 

Seventh: The old jets have faded, and are at most visible in the low radio frequency 
bubbly structures, such as seen for the Virgo cluster region around M87 (Owen et al. 2000). 
The feeding is slowing down, and there is no longer an observable accretion disk, but probably 
only an advection-dominated disk. However, a powerful jet is still there, although below or 
even far below the maximal power. The feeding is still from the residual material stemming 
from the merger. 

Eighth: The feeding of the black hole is down to catching some gas out of a common 
red giant star wind as presumably is happening in our Galactic center. This stage seems to 
exist for all black holes, even at very low levels of activity (e.g., Perez- Fournon & Biermann 
1984, Elvis et al. 1984, Nagar et al. 2000). 

If this concept described here is true, then the superdisk radio galaxies should have large 
outer distortions in their radio images, that may be detectable at very high sensitivity, as 
they should correspond to recently active Z-shaped sources. Also, the superdisk should be 
visible in X-rays, although if the cooling is efficient the temperature may be relatively low. 
Table 2 suggests that the merger is imminent, if the precession of the jet is measurable within 
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a few years, and the opening angle of the precession is much narrower than the wind cone, 
reflecting the earlier longer time precession (see Gopal-Krishna et al. 2007). Therefore, with 
very sensitive radio interferometry it might be possible to detect the underlying jet despite 
its rapid precession, although immediately before the actual merger the feeding of the jet 
will be turned off. 

As more and more pieces of evidence suggest that AGNs are the sources of ultra-high 
energy cosmic rays (Biermann & Strittmatter 1987, Biermann et al. 2007) we need to ask 
what we could learn next. Clearly, after a spin-flip, the new relativistic jet bores through a 
new environment, with lots of gas, and so suffers a strong decelerating shock. In such a shock 
particles are accelerated to maximal energies, and at the same time, as they leave the shock 
region interact with all that interstellar gas. Therefore, such sites are primary sources for any 
new particles, such as high energy neutrinos (Becker et al. 2007). Such discoveries may well 
be possible long before we detect the low frequency gravitational waves from the black hole 
merger. As at such high energy neutrinos travel straight across the universe, and suffer little 
loss other than from the adiabatic expansion of the universe, the black holes resulting from 
a merger of two black holes, with subsequent spin-flip, will be primary targets for searches 
for ultra-high energy neutrinos, and perhaps other photos and particles at extreme energies. 



7. Concluding Remarks 

Whereas it has been questioned in the past whether the central SMBHs of merging 
galaxies will be able to actually merge or their approach will stall (due to the process of 
loss-cone depletion) at a distance where dissipation through gravitational radiation is not 
yet efficient (for a review of these considerations see Merritt & Miloslavljevic 2005), the role 
of the dynamical friction as bringing close the SMBHs to the transition radius, from where 
gravitational radiation undertakes the control of the dissipative process has been recently 
confirmed (Zier 2006) and also complementary mechanisms were proposed (Alexander 2007, 
Sesana et al. 2006, 2007a,b, Hayasaki 2008). The space mission LISA is predicted to 
detect the merger of SMBHs. The statistical arguments of Rottmann (2001), using radio 
observations, suggest that all strong central activity in galaxies may involve a merger of two 
black holes. Therefore, we have assumed in this paper that whenever galaxies merge, so 
will do their central SMBHs. Even if there would be exceptions under this rule, this would 
reflect only in the inclusion of an overall factor < 1 in the number of mergers of SMBHs as 
compared to the number of mergers of galaxies, derived in Section 2, which would not affect 
the mass ratio estimates of our paper. 

Guided by reasonable and simple assumptions we have shown that binary systems of 
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SMBH binaries formed by galaxy mergers typically have a mass ratio range between 1/3 and 
1/30. Following this, we have proven that for the typical mass ranges a combination of the 
SO precession and gravitational radiation driven dissipation produces the spin-flip of the 
dominant black hole already in the inspiral phase, except for the particular configuration of 
the spin perpendicular to the orbital plane. During this process the magnitude of the spin is 
unchanged, therefore the merger of a high spin (and high rotation parameter) black hole with 
the smaller black hole results in a similar high spin state at the end of the inspiral phase. 
These are the main results of our paper. 

There is a related discussion undergoing in the literature, whether the high spin of 
SMBHs is produced by prolonged accretion phases or by frequent mergers. Even a scenario, 
where the SMBHs have typically low spin (King & Pringle 2006) was advanced, based on the 
assumption of short periods of small accretion from random directions. Hughes & Blandford 
(2003), extrapolating results from v = q^ 1 <C 1 binaries to comparable masses, have shown 
that mergers spin-down black holes. Volonteri et al. (2005) have studied the distribution 
of SMBH spins under the combined action of accretion and mergers, and found that the 
dominant spin-up effect is by gas accretion. Recently, Berti & Volonteri (2008) have consid- 
ered the problem of mergers by taking into account improvements in the numerical general 
relativistic methods (Pretorius 2007), and a recent semianalytical formula, which gives the 
final spin in terms of the initial dimensionless spins, mass ratio, and relative angles of orbital 
angular momentum and spins (Rezzolla et al. 2008a,b,c, Barausse & Rezzolla 2009). They 
have found that mergers can result in a high spin end state only if the dominant spin is 
aligned with the orbital angular momentum of the system (thus the smaller mass orbits in 
the equatorial plane of the larger). Their considerations extend from comparable masses to 
mass ratios of 1/10. However, Berti and Volonteri (2008) neglected the angular momentum 
exchange and transport between black hole, jet, and inner accretion disk by magnetic fields 
(see, e.g., Blandford 1976); this may modify or even sharpen the conclusions. 

We can add three remarks to this discussion. First, we have shown by analytical means, 
that for the typical mass ratio range the inspiral phase ends with a considerably lower value 
of the orbital angular momentum compared to the spin (see the last picture in Figured]). A 
heuristic argument then shows that such a small angular momentum could not significantly 
change the direction of the spin during the next phases of the merger. Apart from this 
small orbital angular momentum, the problem being axially symmetric, we do not expect 
significant further spin-flip due to gravitational radiation in the last stages of the inspiral. 

Second, the configuration of orbital angular momentum aligned with the dominant spin 
is not a preferred one in the gravitational radiation dominated post-Newtonian regime. It is 
not clear yet whether such an alignment could be the by-product of previous phases of the 
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inspiral, when dynamical friction (Zier & Biermann 2001), three-body interactions (Sesana 
et al. 2006, 2007a,b), relaxation processes due to cloud-star interactions (Alexander 2007), 
three disk model accretion (Hayasaki 2008), and other possible mechanisms occur. Since the 
stellar system is often slightly flattened, differential dynamical friction could produce the 
near alignment necessary to allow very high spin after a merger. 

Third, the magnitude of the spin is practically unchanged in the inspiral phase, discussed 
here. This is because the loss in the spin vector by gravitational radiation, a second PN order 
effect, calculated from the Burke-Thorne potential (Burke 1971), is perpendicular to the 
spins, yielding another precessional effect (Gergely et al. 1998c). Below ISO this estimate 
should break down, as indicated by numerical simulations reporting on various fractions of 
the spin radiated away. In this context we want to emphasize the unchanged magnitude of 
the spin during the inspiral, as important initial data for the numerical evolution during the 
plunge and ring-down. 

We also mention here the results of the numerical relativity community showing a con- 
siderable recoil of the merged SMBH in particular cases, mostly for equal masses and peculiar 
configurations of the angular momenta (Briigmann 2008, Gonzalez et al. 2007a, b, Koppitz 
et al. 2007). It has also been shown that the recoil regulates the SMBH mass growth, as 
the SMBH wanders through the host galaxy for 10 6 ~ 10 8 years (Blecha and Loeb 2008). 
According to the empirical formula of Campanelli et al. (2007a, see also Lousto & Zlochower 
2009) the recoil velocity scales with q~ 2 / (1 + q~ r ) 2 (1 + q~ r ), which for q~ l — v <C 1 reduces 
to a scaling with q~ 2 . Therefore, we do not expect significant recoil effects in the typical mass 
ratio range of the SMBH mergers. 

We suggest that the precessional phase of the merger of two black holes, occurring 
prior to the spin-flip, is visible as a superdisk in radio galaxies (Gopal-Krishna et al. 2007). 
The precessing jet appears as a superwind separating the two radio lobes in the final stages 
of the merger. According to our model such radio galaxies are candidates for subsequent 
SMBH mergers. Further observations and theoretical work may be capable of identifying 
such candidates likely to merge, and determine the timescale for this to happen. The restart 
of powering a relativistic jet (after the spin-flip and merging) will produce ultra-high energy 
hadrons, neutrinos and other particles. 

Based on the estimates given in Table 2 for the precessional and inspiral timescales, 
we can say the following. If we were to observe a precession timescale of three years in a 
superdisk radio galaxy, we could confidently predict a plunge in about 300 years, which should 
be observable. Faster precession timescales would take some effort to identify. However, if 
we were able to even identify a precession timescale of days to weeks, then the plunge would 
be predicted to happen a few months to a few years thence: powerful gravitational waves at 
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very low frequency would then be emitted. 

The picture developed here differs from that in Wilson & Colbert (1995) in that we do 
not identify just the rare mergers of two massive black holes of about equal masses with 
radio galaxies and radio quasars. We intend to revisit the interactions with the stars (Zier 
et al. 2009, in preparation), discuss the spin of the black holes in another work (Kovacs 
et al. 2009, in preparation) developed from Dut-an Sz Biermann (2005), finally to work out 
quantitatively the relation of the merger of black holes and the statistics of radio galaxies 
(Gopal-Krishna et al. 2009, in preparation). 
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